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Preheating and other particle production phenomena in the early Universe can give rise to high- 
energy out-of-equilibrium fermions with an anisotropic stress. We develop a formalism to calculate 
the spectrum of gravitational waves due to fermions, and apply it to a variety of scenarios after 
inflation. We pay particular attention to regularization issues. We show that fermion production 
sources a stochastic background of gravitational waves with a significant amplitude, but we find 
that typical frequencies of this new background are not within the presently accessible direct detec- 
tion range. However, small-coupling scenarios might still produce a signal observable by planned 
detectors, and thus open a new window into the physics of the very early Universe. 
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Introduction. The existence of gravitational waves 
(GW) is arguably one of the most important predic- 
tions of the general theory of relativity that still re- 
mains unverified. The measured decay of the orbital pe- 
riod of compact binaries JJ has nevertheless provided a 
strong, though indirect, evidence for GW. On theoreti- 
cal grounds, we expect that the universe today should 
be permeated by a variety of GW backgrounds of diverse 
origin. For instance, GW are expected from astrophysi- 
cal sources, like the collapse of supernovas or the coales- 
cence of compact binaries. They are also expected from 
high-energy phenomena in the early Universe, see [2] for 
a review. Many plans for direct detection experiments 
exist, such as the VIRGO interferometer, the Laser Inter- 
ferometer Gravitational Wave Observatory (LIGO), the 
European Laser Interferometer Space Antenna (eLISA), 
the Big Bang Observer (BBO), or the Decihertz Inter- 
ferometer Gravitational Wave Observatory (DECIGO); 
all these observatories operate at some typical frequen- 
cies ranging from 10^^ Hz to 10'* Hz, and it is highly 
expected in the community that GW should be detected 
within this decade. 

A number of constraints on GW have been derived 
from a variety of considerations related to things such 
as millisecond pulsars [Hj and Big Bang nucleosynthe- 
sis [4]. Indeed, there are hopes that cosmology may in 
near future offer more insight on the existence and na- 
ture of GW. From the observations of the Cosmic Mi- 
crowave Background temperature and polarization [5] 
anisotropics, one can already infer upper bounds on the 
amplitudes of GW (see also [6]). 

During inflation metric perturbations, including tensor 
modes, i.e. GW, are generated. If the scale of inflation 
is sufficiently high, the GW background from inflation 
could be detected directly [7] with satellite missions such 
as the Planck Surveyor or the proposed CMBpol satel- 
lite. Moreover, post-inflationary dynamics can also be a 
source of stochastic backgrounds of GW, generated by 
causal mechanisms very different from the quantum na- 
ture of the inflationary GW. After inflation the energy 
stored in the scalar inflaton field responsible for the su- 
perluminal expansion is converted into (almost) all the 



matter and radiation of the universe. This stage is called 
reheating, by the end of which the inflaton decay prod- 
ucts have thermalized among themselves and the stan- 
dard hot Big Bang evolution can commence. The pro- 
cess of reheating is not well understood. However, it is 
often assumed to be initially driven by non-perturbative 
effects, consisting in violent bursts of particle production 
known as preheating [8,. 

In this letter we will focus on the generation of GW 
right after inflation during preheating or other stages 
of non-perturbative particle production. Most previ- 
ous studies of post-inflationary phenomena generating 
GW have focused only on bosonic sources [5HTT]. usu- 
ally scalar fields. Here we want to complete the picture 
by considering fermionic fields as the source of GW. The 
phenomena in which fermions are created after inflation 
correspond to out-of-equilibrium periods in the evolution 
of the universe. Thus the created fermions have typically 
a non-thermal spectra and, as a consequence, they have 
a non-trivial anisotropic stress which will source GW. 
We will consider the possibility that the inflaton couples 
to fermions, and discuss the production of GW during 
fermionic preheating [T31[T3]. Our considerations also ap- 
ply to situations where fermions are produced by a scalar 
held other than the inflaton. The fermions we observe in 
Nature must have been generated sometime between the 
end of inflation and Big Bang nucleosynthesis. Fermionic 
preheating, or production of fermions in general, is thus 
no less natural than the usual bosonic production, but it 
is technically much more difficult to treat. 

Since sources of GW produce spectra with distinctive 
shapes and amplitudes, it is important to characterize all 
the potential sources. The question then is: how big is 
the amplitude, and what is the frequency, of the stochas- 
tic GW produced by fermions in the early universe? The 
aim of the present letter is to answer these questions. 
We first develop a general formalism for computing the 
GW spectrum generated by an ensemble of fermions, 
and then apply it to two distinct scenarios: 1) fermions 
generated in preheating, and 2) thermal Universe into 
which the fermions are injected. From now on we will 



work in units h 



1, with M„ w 2.4 X 10^^ GeV the 



reduced Planck mass. Summation will be assumed over 
repeated indices. 

Gravitational waves from fermions: the formal- 
ism. Let us consider a flat Friedman-Robertson- Walker 
background with a{t) the scale factor and t conformal 
time. GW are the transverse-traceless (TT) part of met- 
ric perturbations hij, 

ds^ = a^{t){-dt^ + ((5y + h,j)dx'dx^), (1) 

subject to dihij — ha ~ 0. The spectrum of energy 
density of a stochastic GW background in comoving mo- 
mentum k is given as 

"" (M) ^ ^fS#. (^) 



dlogfc^"'"' Sir^a^it) ' 

where |/ifc(i)p is the power spectrum of hij — -^. For 
initial conditions hij{ti) = hij{ti) = 0, the GW spectrum 
at sub-horizon scales becomes [2] 

^^^^{k,t) = f , ,\ , f dt' f dt"a^{t')a\t") 

X cos{k{t'~t"))n^k,t',t"), (3) 

where 11^ is the unequal time correlator (UTC) of the 
TT-part of the anisotropic-stress 11, 



TT 



{IlJ^{k,t)Ilj;^\k',t')) = i2nfn^ik,t,t')Snik-k') (4) 

Let us now consider spin-^ fermions as the source of 
GW, which can be represented as 



dk 



ikx 



|ak,rUk,r(i) + 6Lk,r'^k,r(i)| , (5) 



with r = 1, 2, 5i = -5_2 = (1 0)T, 52 = S-i = (0 1)^, 
and Ur^br the usual creation/annihilation operators 
obeying the relations {ar(k), aj(q)} = {&,.(k), &J(q)} ~ 
i2TTfSrsSD0i " q), {ar{k),bl{(i)} = 0. The fermion 
energy-momentum tensor is given by J15j 



i;. 



(x,0 = ^^ (V'7(*^j)^ - '0£'(»7j)V') 



(6) 



with Di = di + \['ya,'yp]<^'j the covariant derivative, ji 
the standard (flat-space) Dirac matrices, and uj"^ the 
spin connection. The TT part of the anisotropic-stress 
in Fourier space, can be simply obtained by means of the 



orthogonal projector Vi 



kikj, as Uj^iKt) 



ViiTimTrnj " ^VijVimTim- We then have all the ingredi- 
ents to calculate n^(fc, t, t') and obtain the GW spectrum. 
Using the fact that Vp.± ~ ±u* -p, we find 



dlogfc^ ' ^ " 8^4 



^Jdpd9p^sm'0Fik,p,e), 
F{k,p,e) = |/|^)-/i^)|2 + |/|^)-/i^)|2 (7) 



dt' 



I^±\k,p,e) = I cos{kt')uu-p.±{t'W±{t'), 



with I^ defined analogously by replacing cos(fci) by 
sin(fci). Eq. (l7| is the master set of formulae that describe 
the GW spectrum at subhorizon scales, as generated by 
some fermionic field ijj with eigenfunctions Uk.±(0. The 
eom for Wk,±(i) will follow from the Dirac equation. For 
any process in the early Universe where fermions are ex- 
cited, one just needs to plug in the solutions Uk,± (t) into 
the master equation (ItI) to find the spectrum of GW. 

The structure of the formulae in Eq. ([7| resembles that 
of scalar fields sourcing GW. However, in the bosonic 
case, apart from multiplicative factors, there appears a 
power p^ instead of p^ in the integrand, there are no po- 
larization indices -t-, — , and of course the fermionic mode 
functions Wk,±(i) are replaced by the Klein-Gordon scalar 
modes (j)k{t). 

Both bosonic and fermionic vacuum expectation val- 
ues (VEVs), like 11^, require regularization. In the 
case of bosons, this has not been an issue in the liter- 
ature, since the bosonic UTC are either introduced as 
a theoretical regularized Ansatz, or in the case of lat- 
tice simulations, the ultraviolet (UV) modes causing the 
divergence are simply not captured. In the fermionic 
case one cannot skip regularization. To regularize 11^, 
note first that the VEV of the source itself, Hji^ , needs 
also to be regularized. Similarly to the flat-space case, 
regularization of the source's VEV amounts to a sub- 
straction of the zero-point fluctuations or, equivalently, 
to a time-dependent normal-ordering (tNO) procedure, 
i.e. (n^-^),eg = (0|n^-^|0) - muf^lOt), with |0) the ini- 
tial vacuum and |0t) the vacuum at time t. This removes 
the unphysical divergence in the VEV of 11^-^ at every 
time t. In practice, tNO amounts to the replacement 
(Mp,±Up',±) -> (wp,±Up',±),eg = l/3pll/3p'|wp,±Wp',± - 
(^p/3p/Mp,zpUp',^)* with l^pl = a/1 - IttpP, where a, /3 
are the canonical Bogoliubov coefficients [TS] connecting 
the initial creation and annhilation operators with those 
at time t. To render 11^ finite we just need to replace in 
Eq. I?! the functions I^^ by 



Ac) _ 



dt' 

^COS(H') (uk_p.±(t'X±(t')),eg' (8) 



(s) 

and similarly for X!j_ . In this way the convergence of 
the integration over p is ensured by the suppression of 
the UV divergent modes. 

Stochastic background of gravitational waves 
from fermions produced in the early Universe. 

Several scenarios of the early Universe may create high- 
energy out-of-equilibrium fermions by non-perturbative 
effects, for instance a homogeneous scalar field oscillat- 
ing around the minimum of a potential. In such sit- 
uation, if fermions are coupled to that field, a non- 
perturbative population of their modes (respecting Pauli- 
blocking) takes place [TSJ E]. This generates a non- 
trivial anisotropic-stress that sources a stochastic back- 



ground of GW. In this letter we consider two different 
scenarios: 

I) (p)reheating after inflation: Let us assume an infla- 
ton oscillating around the minimum of its potential after 
the end of inflation. This is the case of chaotic infla- 
tion models with polynomial potentials like F ex 0^ or 
V oc (f)^. In general, the shape of the inflaton potential 
during inflation is however irrelevant for our purposes. 
We will just assume that inflation took place at some 
energy scale -E/, and that afterwards the inflaton oscil- 
lates with a certain frequency and decreasing amplitude. 
The expansion of the Universe will be dictated during 
(p) Reheating by the inflaton energy density. For a poly- 
nomial potential V{(j)) it is well known that the scale 
factor behaves as a power law in time. In the simplest 
case V{4)) — ^uJq<j>'^, a{t) ex t^, so the expansion of the 
Universe is matter-dominated (MD). 

II) Oscillating scalar field in a thermal era: It is also 
possible to consider a scalar field, other than the infla- 
ton, which oscillates coherently around a potential mini- 
mum when the Universe has already entered into a ther- 
mal era. That is the case, for example, of the curvaton 
scenario [I2i, which is an alternative to single-field in- 
fiationary models. After the Universe has reheated, the 
curvaton oscillates with decreasing amplitude and fixed 
frequency. The expansion of the Universe is radiation- 
dominated (RD), driven by the thermal relativistic bath 
of particles with energy density pth oc l/a**, and a{t) ex t. 
We assume oscillations begin at some energy scale Ei 
way above the electroweak scale. 

In both scenarios the oscillatory field behaves as 
a damped oscillator, (f>{t) — ^{t)F{t), with decreasing 
amplitude ^(t) and periodic behaviour F{t + 2-11 /ujq) = 
F{t) < 1, as determined by the choosen potential. We 
will assume that is coupled to some fermion species tp 
via a Yukawa interaction iKJjipip, with h the interaction 
strength. The Dirac equation yields the equation for the 
mode functions Uk ± as 



Wk,± 



+ fk^ + h^a^if ± ih{a(l) + a^)\ 



Uk,± = 



(9) 



We have solved numerically the mode equations ^ for 
the scenarios I and II, scanning for the parameters h, wq 
and El and assuming 1^(0) = I'^o'^^ ^nd initial condi- 
tions corresponding to a vanishing particle density. In 
both scenarios, each time passes through the minimum 
of the potential, fermions are created out-of-equilibrium 
by non-perturbative parametric effects [8 . 

Pauli blocking prevents the fermion occupation num- 
bers to grow arbitrarily. Fermion excitations are forced to 
fill up a "Fermi-sphere" in momentum space of comoving 
radius kp ~ a{t)<l^^'^<^o^ with q = /i^($/a;o)^ the resonant 
parameter [13j . The occupation numbers with momenta 
smaller than kp oscillates continuously between and 1, 
whilst the radius kp grows as ex a^'**. When the energy 
density of the created fermions has grown to a significant 



fraction of the energy of the oscillating field - usu- 
ally within tens of oscillations -, the creation of fermions 
finally ceases. The excitations of the fermionic modes 
within the Fermi-sphere source the generation of GW. 
The non-excited modes, i.e. those outside the sphere, are 
on the contrary responsible for the UV divergence dis- 
cussed before. The regularization scheme we impose fil- 
ters the infrared (IR) modes within the sphere k < kp 
and removes the contribution from the UV ones. 

It is useful to define the times i/, i* and iRDi as the 
initial time, the end of GW production, and the first mo- 
ment when the Universe becomes RD, respectively. In 
the scenario I, i/ < i* < iRD, so in this case the ef- 
fective equation of state p/p = w between tj and ^rd 
is typically different from that of RD, wrd = 1/3 (un- 
less V{(p) ex (j)*). It is then convenient to introduce the 
factor e = {a^,/a]xu)'^•^~^^K see next. In the scenario II, 
tRD < ti < t^., so that e — 1. Below the Planck scale, 
GW decouple immediately after production, so we can 
evaluate the GW energy density spectrum today from 
the spectrum computed at the time of production. We 
just have to redshift the amplitude and wavenumbers, 



1/4/ a/ 




X 5 • 10^° Hz 



h'n^..., = 



Pc\d\ogkJ^ 



/}2q 



p\d\ogkJ^ 



i,ad ~ 4xl0-^ p* ^ Ej{ai/a^f^^+'^'^ the energy 
density at i*, and {go/g*Y^^ ^ 0(0.1) the ratio of the 
number of relativistic dof today to those active at i* . 

In Fig. [1] we show two examples of GW energy den- 
sity spectra, obtained with the machinery presented in 
the previous section. The spectra grows in the IR region, 
reaches a maximum at fc ~ g^'^wo, and finally falls down 
in the UV region. The spectra are peaked, as expected 
on dimensional grounds, around the scale ~ kp, char- 
acteristic of the fermions' dynamics. From the fc ^> 
limit in Eq. ([7|, we would expect the spectra to grow as 
ex fc^ in the very IR region. Instead we observe a be- 
haviour ~ /c^, probably signaling that we are exploring 
fc's too close to the maximum, and that one should look 
for much smaller fc's to probe the k^ behaviour. Un- 
fortunately, due to computer limitations, both the limits 
k ^ kp and k ^ kp are indeed very challenging to probe, 
so we leave a complete characterization of the IR and UV 
tails of the GW spectra for a future publication. 

With our formalism we can obtain nevertheless the 
most important aspects characterizing any GW spec- 
trum: the spectral shape around the peak, and the am- 
plitude and position of such peak. For instance in the 
scenario I with Ej - lO^^ GcV and luq - 10^'' GeV, we 
find for q = 10^, today's peak frequency and amplitude 
as / ~ 10^ Hz, h^riaw '^ 10-^°, and for q = 10^ as 
/ - 10^° Hz, h^naw - 10~". In the scenario II, for the 



1 ■iPs» 
p+ dlogA- 



FIG. 1: Spectra of GW at the time of production, in the 
scenario I with Ei ~ lO^** GeV and wq ~ lO" GeV, for q = 
lO'^ (dashed line, peaked at A; = Sq^'^wo) and q = 10® (solid 
line, peaked at fc = Gq^''^ujo). The spectra in the scenario II 
for the same parameters look qualitatively the same. 



same Ej and wq, we find / ^ 10^ Hz, h^Vlcw ~ 10 ^^ 
for q = 10^, and / ~ lO^o, K^naw ~ IQ-^^ ^j. ^ ^ iq6^ 
Thus the amplitude of the GW background explored here 
can indeed be very significant. However the typical fre- 
quencies are too large as compared to the range probed 
by GW observatories (~ 10"^ - 10^ Hz). From Eq. (Il0| 



1/4 



we learn that if wq/P/ ^ 1, we could decrease sig- 
nificantly the frequency towards the observable window. 
However 51gw is suppressed by (ujq/pi), so the ampli- 
tude would be far too small, flaw ^ 10~^°. We hope to 
return to a systematic exploration of the parameters in 
a future publication. 

There are some indications for scenarios where a 
GW background at sufficiently low frequencies and high 
enough amplitudes could be found. For instance, in hy- 
brid inflation, the frequency ujq should be replaced by 
VXv, and the initial energy scale Ej by ^ X^^'^v, where 
A and v are the self-coupling and VEV of an auxiliary 
field coupled to the infiaton. The frequencies would then 
scale as / oc A^/* [TQ so that for sufficiently small A, 
one could obtain a peak frequency within the observable 
range. This possibility can be explored by using the ma- 
chinery developed in this Letter. 

Summarizing, we have shown that fcrmions in the early 
Universe may be very efficient generators of GW. These 
waves remain decoupled since the moment of their pro- 
duction, and thus the amplitude and shape of their spec- 
trum probes the physics responsible for their generation. 
The characteristic spectrum is different from other back- 
grounds of GW, like those arising from binaries coalesc- 
ing [5], which are decreasing with frequency, or those 
arising from inflation or self-ordering flelds |16j . which 
are flat. A comparison of the exact shape of the spec- 
tra predicted here versus those generated by scalar fields 
in phase transitions and (p)reheating requires further 
study. Here we just want to emphasize that, contrary to 
naive expectations based on the Pauli principle, we found 



that fermions are capable of producing a stochastic back- 
ground of GW with a very large amplitude. Although it 
may be manifest only at very high frequencies much be- 
yond present/planned GW detectors' sensitivities, there 
is however some hope that small-coupling models could 
give rise to a signal in observable frequency ranges. 
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